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ABSTRACT 

Recent studies of the molecular medium in nearby galaxies have provided mounting evidence that 
the molecular gas can exist in two phases: one that is clumpy and organized as molecular clouds 
and another one that is more diffuse. This last component has a higher velocity dispersion than the 
clumpy one. In order to investigate these two molecular components further, we compare the fluxes 
and line widths of CO in NGC 4736 and NGC 5055, two nearby spiral galaxies for which high - quality 
interferometric as well as single-dish data sets are available. Our analysis leads to two main results: 1) 
Employing three different methods, we determine the flux recovery of the interferometer as compared 
to the single-dish to be within a range of 35-74% for NGC 4736 and 81-92% for NGC 5055, and 2) 
when focusing on high (SNR>5) lines of sight, the single-dish line widths are larger by ~ (40 ± 20)% 
than the ones derived from interferometric data; which is in agreement with stacking all lines of sight. 

These results point to a molecular gas component that is distributed over spatial scales larger than 
30"(~lkpc), and is therefore filtered out by the interferometer. The available observations do not 
allow us to distinguish between a truly diffuse gas morphology and a uniform distribution of small 
clouds that are separated by less than the synthesized beam size (~3" or ~ 100pc), as they would 
both be invisible for the interferometer. This high velocity dispersion component has a dispersion 
similar to what is found in the atomic medium, as traced through observations of the HI line. 

Keywords: galaxies: ISM — ISM: molecules — ISM: clouds - radio lines: galaxies 


1. INTRODUCTION 


In the classical picture of the interstellar medium 
(ISM), molecular gas is thought to be organized mostly 
inside giant molecular clouds (GMCs), typically traced 
by CO emission, with temperatures of 10 — 20 K and 
int ernal velocit y disp ersions of ~ 2 — 8 km s _1 (e.g., 


Bolatto et al. 20081. This simple picture, however, 
does not hold at the Galactic level, were different stud¬ 
ies have shown that diffuse CO emiss i on is pervasive 


within our Galaxy (Polk et al. 1988 Liszt & Lucas 
1998 Goldsmith et al.' 2008 and references therein). 
Based on studies of molecular gas in our Galaxy, molec¬ 
ular clouds are generally classified as either diffuse 
(A y ~ 0.2, 30 < T(K) < 100), translucent (Ay ~ 1 — 2, 
15 < T(K) < 50), or dense ( Ay ~5-10, 10<T(K)<50) 
(e.g., Snow fc McCall|[2006 ). On the one hand, studying 
the ISM in the Galaxy might seem advantageous from the 
sensitivity and resolution points of view. Galactic studies 
can carry out observations of many different molecules, 
isotopologues, and transitions; while extragalactic stud¬ 
ies typically detect only the most abundant molecule in 
its brightest transition. On the other hand, due to pro¬ 
jection effects and uncertainties in distance determina¬ 
tions, it is difficult to quantify the amount of gas present 
in this diffuse molecular component within our Galaxy. 
Nearby galaxies offer the opportunity to study this sug- 
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gestive additional molecular gas component from an out¬ 
side perspective. 

Molecular gas is preferentially observed through CO 
emission, which can be mapped both with single-dish 
telescopes and with interferometers. Single-dish tele¬ 
scopes, especially after the development of multi-beam 
heterodyne receiver arrays, provide a large instantaneous 
field of view and good sensitivity, although at limited 
spatial resolution. Interferometers, on the other hand, 
provide higher spatial resolution, with the drawback that 
not all spatial scales can be recovered. The largest spatial 
scale to which an interferometer is sensitive to, is deter¬ 
mined by the shortest baselines (i.e., the closest projected 
distance between two individual telescopes in the plane 
perpendicular to the line of sight). Thus, any emission 
arising from structures larger than this spatial scale is 
‘invisible’ to the interferometer. If most gas inside galax¬ 
ies were distributed in structures (e.g., molecular clouds) 
that are smaller than this largest spatial scale (and that 
are separated from each other by at least the synthe¬ 
sized beam size), then either instrument would detect the 
same gas component, only with different spatial resolu¬ 
tion. Therefore, the flux recovered by the interferometer 
depends on the gas distribution and the molecular cloud 
spatial scales present in the observed galaxy as compared 
to the specific resolution of the interferometer’s configu¬ 
rational 


Pety et al. (2013) conducted a study of the molec¬ 
ular gas phase as traced by the CO (1-0) emission in 
M51. They used single-dish IRAM 30m data, Plateau 
de Bure Interferometer (PdBI) data, and a combination 


5 During the image deconvolution the flu x recove r y will also 
strongly depend on the noise properties (e.g., |Helfer et al.|2002j). 
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of both data sets, to compare flux recovery and veloc¬ 
ity dispersions measured in these three data sets. They 
show that the interferometer recovers only (50 ±10)% 
of the total flux, indicative of a molecular component 
that is missed by the interferometer. In addition, they 
found that the purely single-dish velocity dispersions 
are typically twice as large as the interferometric veloc¬ 
ity dispersions. As stated before, molecular gas might be 
present in different cloud types, with different tempera¬ 
tures, densities, and spatial distribution. These different 
physical configurations would modify the expected ve¬ 
locity dispersions present in the gas. Thus, their result 
implies not only that ~ 50% of the gas is missed by the 
interferometer, but that this gas gives rise to larger ve¬ 
locity dispersions compared to the velocity dispersions 
measured for the molecular gas (i.e., CO emission) inside 
the GMC sample (which the interferometer is sensitive 
to). 

In a previous paper, [Caldu-Primo et ah (2013]) inves¬ 
tigated the differences in velocity d ispersion of atomic 
(Hi) gas (from the THINGS survey, |Walter et al.||2008l 
and of molecular gas (from the HERACLES CO survey, 
Leroy et al.|2009) in a sample of 12 nearby spiral galaxies. 
On spatial scales of ~ 1 kpc we measured similar velocity 
dispersions for both HI and CO (as measured by a sin¬ 
gle - dish telescope), and interpreted this as the existence 
of a high - dispersion molecular gas component. 

Based on these two previous results, we deem impor¬ 
tant to rigorously investigate the presence of the molecu¬ 
lar gas component giving rise to these large velocity dis¬ 
persions and determine its mass fraction as compared to 
the molecular gas inside molecular clouds. The existence 
of such a molecular gas component would require that 
CO single-dish and interferometric observations be an¬ 
alyzed and interpreted using different assumptions about 
the origin of the emission. For example, studies of the re¬ 
lationship of molecular gas and recent star formation rate 


in both nearby or d istant galaxies (e.g., review by Kenni- 
cutt & Evans| 2012) implicitly assume that molecular gas 
as traced by CO observations is inside molecular clouds 
and participates in the current star formation process; 
whereas the assumption of a significant diffuse molecular 
component would allow for some of the observed molecu¬ 
lar gas being not directly (or at all) related to the current 
star formation activity. In this study we will further in¬ 
vestigate the presence of a wide-spread, high-velocity 
dispersion molecular gas component by comparing CO 
line widths as measured from interferometric and sin¬ 
gle - dish observations, as well as the fluxes measured by 
means of both instruments. Unfortunately, high sensitiv¬ 
ity interferometric maps of nearby galaxies are still very 
scarce, a situation that will however soon change with 
the advent of ALMA. Thus, we concentrate our current 
study on two galaxies: NGC4736 and NGC5055. We 
comp are single d ish C O data from the HERACLES Sur¬ 
vey (Leroy et al.||2009) a nd the “Nobeyama CO Atlas of 
Nearby Spiral Galaxies” (|Ku no et al.||2007|) to interfero¬ 
metric observations taken by|La V igne| (|2010j) using the 
Combined Array for Research m Miliimeterwave Astron¬ 
omy (CARMA). 

The paper is organized as follows: In Section [2] we 
discuss the main characteristics of the employed data sets 
and the two galaxies studied. In Section [3] we discuss our 
methodology, to finally show the results in Section |4j A 


discussion and summary is provided in Section [5j 
2. DATA AND SAMPLE 

In order to carry out the analysis, the data need to 
meet three conditions: high sensitivity, a clear contrast 
between the scales recovered by the interferometer and by 
the single-dish, and available Hi velocity maps. Ev en 
though there is a large survey, BIMA SONG (Heifer 
et al. 2003), in which the CO (1-0) molecular emission 
in nearby galaxies has been mapped with an interferom¬ 
eter, its sensitivity and resolution render it not suitable 
for this specific study. Therefore, we limit our study to 
the 2 galaxies for which high - sensitivity and resolution 
CARMA interferometric maps exist. 

2.1. CO(l 0) Data 

2.1.1. Interferometric Data 

The interferometric CO data come from the CARMA 
interferometer. Both galaxies are mapped in the 
CO (1 0) line covering the inner ~ 3' diameter of the 
targets via a 19-point mosaic. Observations were taken 
between January 2007 and July 2008 using the C, D, 
and E configurations resulting in a synthesized beam size 
of 3.51 x 2.97" for NGC4736 and of 3.24 x 2.81" for 
NGC 5055. The correlator was setup with two overlap¬ 
ping 62 MHz spectral windows to yield a spectral resolu¬ 
tion of 1.95 MHz (or 5.08 km s -1 at 115 GHzf*] 

The CARMA data were calibrated usingMIRIAD 
(|Sault et aLl[l995|), weighted by noise variance, and ap¬ 
plied robust weighting. The data were cleaned down to 
a cut-off of 1.5 times the theoretical rms noise in each 
pixel using MOSSDIA which per forms a steer CLEAN on 
a mosaicked image (|Steer et al.||l984|). The images are 
primary - beam corrected. The flux scale was determined 
by observatio ns of flux standards, including Uranus and 


Mars. In La Vigne (2010) a full description of the data 
reduction can be found. In their analysis they compare 
the CARMA observations with the BIMA SONG data 
(Heifer et al. 2003). After verifying the flux scales of 


both datasets they concluded that no flux rescaling was 
required. The rms noise sensitivity at 10 km s” 1 reso¬ 
lution is 27 mjy beam -1 (237.5 mK|^ for NGC 4736 and 
19mJybeam _1 (193.2mK) for NGC 5055. 

2.1.2. Single - dish Data 

We also use single-dish CO (1 0) mapping from_the 
Nob eyama CO Atlas of Nearby Spiral Galaxies” (Kuno 
et al. 2007), which encompasses 40 nearby spiral galax- 
les at distances smaller than 25Mpc. The observations 
were carried out with the Nobeyama 45m telescope and 
typically cover most of the optical disk with a spatial 
resolution of 15". 

In the case of NGC 4736, the recei ver used was BE ARS, 
which consists of 25 beams (Sunada et al. 2000). In 
this case digital spectrometers were used as backends. 
The total bandwidth and frequency resolution were 
512 MHz and 605 kHz, which at 115 GHz correspond to 
1331km s -1 and 1.57km s -1 , respectively. 

6 The now publicly available data for NGC 4736 has 10 km s -1 
since the original data were unfortunately lost, and cannot be re¬ 
covered. 

7 All sensitivities in K are stated using the beam size of the 
original data. 
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In the case of NGC5055, a 4-beam SIS receiver 
was used, together with acousto - optical spectrometers 
(AOS) used as backends. This configuration yields a 
frequency resolution and channel spacing of 250 kHz 
and 125 kHz, respectively, providing a total bandwidth 
of 250 MHz. At 115 GHz this corresponds to veloc¬ 
ity resolution and velocity coverage of 0.65 km s -1 and 
650 km s -1 , respectively. 

In both cases t he s pectra were then smoothed to 
5 km s -1 (Nakai et al. 1994) as part of the data re¬ 
duction. This Is" the spectral resolution of the data 
we use. At 10 km s _1 resolution, the rms noise sen¬ 
sitivity is 69mJybeam _1 (28.3 mK) for NGC4736 and 
172 mJy beam " 1 (70.7 rnK) for NGC 5055. 


2 . 2 . CO (2-1) Single-dish Data 

W e use CO (2-1) data from the HERACLES sur¬ 
vey (Leroy et al. 2009). This survey used the IRAM 
30m telescope to map molecular gas from 48 nearby 
galaxies inside the optical radius (^ 5 ). This survey 
used the HERA (a 9-beam dual polarization hetero¬ 
dyne receiver array) together with the Wideband Line 
Multiple Autocorrelator (WILMA) backend. WILMA 
consists of 18 units of 2 MHz channel width each 
and 930 MHz bandwidth, yielding a velocity resolu¬ 
tion of 2 . 6 km s _1 in a 1200 km s -1 velocity range at 
230 GHz. The spatial resolution obtained at this fre¬ 
quency is 13" after gridding. The rms noise sensitivity at 
10 km s -1 is of 90mJybeam _1 (11.5 mK) for NGC 4736 
and 94mJybeam _1 (12 mK) for NGC 5055. 


2.3. Hi Data 

For our analysis we use Hi intensity weighted veloc¬ 
ity jnupsWThesemaps come from the THINGS survey 
(Walter et al. 2008). This survey encompasses 34 nearby 
galaxies at distances of 2—15 Mpc. The observations 
for the two galaxies we analyze have high spectral (5.2 
km s -1 ) and spatial (~ 11 " for natural weighting) reso¬ 
lutions. 


2.4. Sample 

As mentioned in the introduction, our sample is 
limited by the availability of sensitive, wide-held inter¬ 
ferometric data sets. We thus focus on two galaxies 
only, NGC 4736 and NGC 5055. They have the following 
properties: 


NGC4736 is an isola ted spiral galaxy of type SAab 


(J.43 Mpc (Tonry et al. 

i —1 

0 

0 

01 

Blok et al. 

2008), and 

an 0 


(~ 7.5 kpc) (|Kennicutt et al.||2003|). 

NG C 5055 is a tlocculent s piral galaxy classified as 
SAbc (Kennicutt et al. 2003). Its distance is 7.8 ± 


2.3 Mpc ([M asters 2005), its inclination is 59° (de Blok 
et al.||2008|), and its optical radius is R 25 « 6 aremin 


(~ 13kpc) (Kennicutt et al. 


More infomation onTh e ga. 


2003). 


axy properties c an be found 


Kennicutt et al. (2003), Walter et al. (20081, and Leroy 
et al.| ( |2009[ ). We show the integrated intensity CO maps 


of the three different surveys at a common resolution of 
15" in Figure [I] 


3. METHODOLOGY 
3.1. Data Homogenization 

For our analysis we convolve the data to match the 
spatial and spectral resolution in all cases. This common 
resolution is defined by the poorest resolution of our data 
sets: the spectral resolution is set by the CARMA data, 
and the spatial resolution is set by the Nobeyama data. 
After hanning smoothing the data to a common spectral 
resolution (which is discussed below ), we use the REGRID 
task in MIRIAD (Sault et al. 1995]) to put all the data 
sets on a common grid. Finally, we convolve the data 
sets to the 15" Nobeyama limiting spatial resolution. 

The limiting spectral resolution is different for the 
two galaxies. NGC 4736 has a CARMA - limited 
lOkrns -1 spectral resolution, while NGC 5055 has a 
CARMA/Nobeyama limited 5 km s -1 spectral resolu¬ 
tion. Therefore, for NGC 5055 we first carry out the anal¬ 
ysis at a 5 km s -1 spectral resolution and then compare 
the results when using a 10 km s -1 spectral resolution. 
We find that our results do not change. Thus, for sake of 
uniformity, we here present the analysis using a common 
spectral resolution of 10 km s -1 for both galaxies. 

In order to correctly Nyquist sample the data, we con¬ 
struct a hexagonal grid with a 7.5" spacing (half of the 
beam size). 

Due to the attenuation of the primary beam in inter¬ 
ferometric observations, the sensitivity drops towards the 
edges of the map. We construct a sensitivity mask (after 
primary- beam correction of the cube) from the inter¬ 
ferometric noise map and apply it to all data sets. In this 
way we make sure: 1) to use the same individual lines- 
of-sight (LOS) in all data sets, and 2) that we only use 
LOS which have homogeneous noise properties in the in¬ 
terferometric data. To do so, we derive the galactocentric 
distance at which the noise in the CARMA observations 
is increased by 30%, and only use this region for further 
analysis. For clarification, we show this region as a solid 
black circle in Figure[2]for each galaxy on top of the HE¬ 
RACLES IRAM 30M integrated intensity map in gray 
shades. Inside the circle we also show our hexagonal grid 
points (blue x symbols). Finally, we also show the 15" 
radial bins used later in this paper as red ellipses. 

From Figure[2]we note that the sensitive region defined 
by the CARMAprimary beam response is independent of 
the inclination of the galaxy. However, the radial bins are 
ellipses because of the inclination. This results in having 
lower covering factors within the ellipses for the outer¬ 
most radial bins. To quantify this effect, we present in 
Table]!] the number of LOS from the hexagonal grid that 
fall within each radial bin and within our mask (third col¬ 
umn), as compared to the number of LOS that would be 
present without taking into account the sensitivity mask 
(second column). 


3.2. Measuring Line Widths for Individual LOS 

We estimate the spectral line widths using single Gaus¬ 
sian fits. In order to obtain reliable measurements for in¬ 
dividual lines of sight (LOS), we require that they have 
high significance. Therefore, we begin by computing the 
line widths for the LOS which comply with the follow¬ 
ing characteristics: a,) lie within the CARMA sensitive 
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NGC 4736 



NGC 5055 



Figure 1 . CO integrated intensity maps of NGC 4736 (top) and NGC 5055 (bottom) at resolution (the beam size is shown at the lower 
left corner). The first column shows the CARMA CO(l-O) interferometric data, the second column shows the HERACLES IRAM 30m 
CO (2 1) maps, and th e thi rd column shows the Nobeyama 45m CO (1-0) maps. The three maps were constructed using the HERACLES 
mask (details in Section |4.5| ). For both galaxies the intensity range goes from 0 to 38 Kkm s —1 , and the 6 contours show intensity levels 
separated by 6 Kkm s —1 , starting at 6 Kkm s -1 . The black circle shows the region where the CARMA cube has homogeneous noise 
properties, i.e., where the rms is within 30% of the value at the center of the map. 


region (black circle in Fig. [I]), fojliave peak SNR > H 
and c) have integrated intensities of at least 10cr, where 
a = A Vy/nrms and n is the number of measured points 
used to fit the line. For each of these high SNR LOS 
we use the Hi first moment value,i.e., the local inten¬ 
sity - weighted velocity, as a proxy for the velocity where 
we expect to find the corresponding CO line (We take 


the Hi data from the THINGS survey (Walter et al. 


20081). This is based on our earlier findings that the 


HI and CO line-of-sight velocities agree well at the 
resolution we are w orking at (Schruba et al.|2011 Caldu- 
Primo et al.||2013|). We attempt to lit a single Gaussian 
to each spectrum using MPFIT (IDL procedure from Craig 
Markwardt) inside a window of 100 km s _1 centered in 
the Hi first moment value. If the fitted Gaussian is 
broader than 100 km s , we proceed to refit with a dou¬ 
ble - horn profile (a G aussian sc al ed b y a symmetric sec¬ 
ond order polynomial Saintonge (2007)). Double-horn 
profiles are mostly encountered in the inner regions of 
galaxies, where bulk motions or beam smearing could be 
dominant. These will be discarded for our analysis (see 
below). The integrated intensity selection criteria (b) re- 


8 This SNR is computed in the 10 km s 1 spectral resolution and 
15" spatial resolution cubes. 


moves LOS with significant negative structures (‘bowls’) 
in the interferometric datcEl 
Th e errors are computed as discussed in detail in 
Caldu-Primo et al. (|2013): for each LOS we add random 


noise to the original spectrum and redo the line width 
measurement. We repeat this procedure a 1000 times. 
The adopted error in the line width measurement is the 
dispersion between the results of these 1000 iterations. 


3.3. Measuring Line Widths for Stacked Spectra 

The main purpose of the stacking is to increase the 
SNR in our data as well as to include all available LOS 
in our analysis. The stacked spectra are the most mean¬ 
ingful measurements, as they represent the luminosity- 
weighted FWHM including all LOS inside a tilted ring 
(and inside the CARMA sensitive region). In order to 
stack individual LOS coherently, we first remove the cir¬ 
cular motions of the four data sets using the Hi first 
moment value. The THINGS survey has very high sen¬ 
sitivity, resulting in reliable velocity maps. Even in the 
case of NGC 4736, which has low content of HI in the cen- 

!J These negative bowls are expected to be present in interfero¬ 
metric data because of the extrapolation required to overcome the 
incompleteness of the u —v coverage. 
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Figure 2. HERACLES IRAM 30M CO integrated intensity maps of NGC 4736 (left) and NGC 5055 (right) at 15” resolution (grayscale). 
The black solid circle represents the area where the interferometric primary beam sensitivity is within 30% of the value of that at the center 
of the map. The blue points inside the circle show the hexagonal grid points that are used for sampling the data. The red ellipses show 
the 15” radial bins used for stacking. From inside out and in units of 125 they are, for NGC 4736: 0.04, 0.10, 0.17, 0.23, 0.29, 0.35, 0.41, 
and 0.46; and for NGC 5055: 0.03, 0.07, 0.11, 0.15, 0.19, 0.23, 0.27, and 0.30. 


Table 1 

Filling factors within each 15” radial bin 


# grid points # grid points in # of significant grid points, SNR > 5 b 

R25 in radial bin CARMA-sens, region (%) a CARMA (%) HERACLES (%) Nobeyama (%) THINGS (%) 


1 

2 


3 

4 


5 


6 


7 

NGC 4736 

0.04 

10 

10 

(100) 

4 (40) 

9 

(90) 

3 

(30) 

6 

(60) 

0.10 

35 

35 

(100) 

28 (80) 

32 

(91) 

20 

(57) 

32 

(91) 

0.17 

57 

57 

(100) 

44 (77) 

57 

(100) 

40 

(70) 

57 

(100) 

0.23 

73 

73 

(100) 

34 (47) 

73 

(100) 

43 

(59) 

73 

(100) 

0.29 

101 

101 

(100) 

12 (12) 

72 

(71) 

24 

(24) 

101 

(100) 

0.35 

116 

99 

(85) 

0 (0) 

43 

(43) 

0 

(0) 

99 

(100) 

0.41 

142 

64 

(45) 

1 (2) 

25 

(39) 

0 

(0) 

61 

(95) 

0.46 

165 

10 

(6) 

0 (0) 

3 

(30) 

0 

(0) 

8 

(80) 


NGC 5055 


0.03 

9 

9 

(100) 

6 

(67) 

9 

(100) 

4 

(44) 

9 

(100) 

0.07 

23 

23 

(100) 

7 

(30) 

20 

(87) 

16 

(70) 

20 

(87) 

0.11 

34 

34 

(100) 

4 

(12) 

34 

(100) 

33 

(97) 

34 

(100) 

0.15 

54 

54 

(100) 

12 

(22) 

54 

(100) 

47 

(87) 

54 

(100) 

0.19 

66 

54 

(82) 

11 

(20) 

54 

(100) 

48 

(89) 

54 

(100) 

0.23 

86 

44 

(51) 

9 

(20) 

39 

(89) 

27 

(61) 

39 

(89) 

0.27 

89 

27 

(30) 

5 

(19) 

27 

(100) 

18 

(67) 

27 

(100) 

0.30 

117 

26 

(22) 

11 

(42) 

24 

(92) 

14 

(54) 

26 

(100) 


a Number of grid points within the CARMA sensitive region in a particular bin. In parenthesis the is the % of these grid points compared 
to the total number of grid points within each radial bin (i.e. compared to Column 2). 

b For each instrument and for each radial bin we show the number of LOS with SNR> 5 and within the CARMA - sensitive region (black 
circle in Fig.Jljl. In parenthesis is the % of these points compared to Column 3. 


tral parts, the peak SNR is larger than 5er in 97.3% of the 
LOS inside 0.5 R^s- Once all spectra are centered in the 
same velocity frame (in this case zero velocity) they can 
be stacked coherently. The stacking procedure, together 
with the error determination, are described in detail in 


Caldu-Primo et al. (2013). After the stacking is done, we 
compute 1the line widths as discussed above. The total 
error in the determination of line widths has contribu¬ 
tions both due to noise (as modeled for the individual 
LOS above) as well as due to the specific selection of 


spectra that have been s tacked (which we assess t hrough 
bootstrapping following Caldu-Primo et ah||2013 ). 


4. RESULTS 


4.1. FWHM as function of Galactocentric Distance 

We can select which LOS we want to stack. Since 
many galactic properties depend on the distance to the 
center of the galaxy, we choose to stack the individ¬ 
ual LOS inside tilted rings of 15" width within each 
galaxy. As a first step, we stack only individual LOS 
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with SNR > 5. We then compare the measurements ob¬ 
tained from these (high SNR) stacked spectra, to the 
median values measured from the individual high SNR 
LOS. We do this comparison to test the validity of the 
stacking method. In Figure[3] we plot the FWHM mea¬ 
surements as function of galactocentric distance for each 
of the 4 different data sets (we include Hi for sake of 
completeness). The FWHM measured when fitting high 
significance (SNR > 5) individual LOS are shown in small 
grey symbols. The FWHM measured when only stacking 
the grey points (i.e., individual LOS with SNR> 5) are 
shown in medium-sized blue symbols, whereas the me¬ 
dian FWHM of the individual high SNR LOS are shown 
in black small symbols (The red points are explained be¬ 
low). We see that at small radii (R< 0.2 R 25 ) we observe 
a rapid increase in the measured FWHM, however, here 
the observed s pectra are si gn ificant l y bro adened by beam 
smearing (see |Caldu-Primo et al.| 2013], for a thorough 
discussion on this effect). The shaded regions highlight 
the radial range where beam smeari ng bro a dens the in¬ 
trinsic line profile by more than 30%. Caldu-Primo et al. 
(2013) show that for a galaxy with characteristics simi¬ 
lar to NGC4736 beam smearing accounts for >30% of 
the measured line width at R<0.15R25 and for more 
than 10% out to 0.25R 25 . In NGC5055, beam smear¬ 
ing accounts for > 30% of the observed line widths at 
R< 0 . 2 R 2 5 , and > 10% at 0.4R 25 . At large radii out to 
0.5 R 25 , where line profiles are only mildly or not at all 
affected by beam smearing, we find roughly constant CO 
velocity dispe rsions, as have been found already in earlier 
studi es (e.g., Tamburro et al. 2009 Caldu-Primo et al. 
20131. In the regions of interest (where beam - smearing 
is not dominant), the measurements obtained for LOS 
with SNR > 5 are in agreement. 

In Table]]] (columns 4-7) we show the number of 
LOS with SNR larger than 5 within each 15" radial bin 
and within the CARMA - sensitive region for each of the 
datasets (in parenthesis we show which percentage of the 
total number of LOS within the CARMA - sensitive re¬ 
gion these high SNR LOS represent). It is important to 
note that we select the individual LOS with high SNR 
independently for each data set, i.e., a high SNR LOS in 
the HERACLES cube does not necessarily have such a 
high SNR LOS counterpart in the CARMA cub^Jsee 
Table]!]). It is clear that by taking only high SNR mea¬ 
surements we are excluding a large percentage of LOS 
from the analysis. 

Therefore, we then proceed to stack all LOS inside the 
15" wide tilted rings. The resulting stacked spectra are 
shown in Figure]!] for NGC4736 and NGC5055, respec¬ 
tively. In these figures we show the normalized stacked 
spectra with peak SNR higher than 5 in all data sets. 
The corresponding FWHM obtained for these stacked 
spectra are plotted in Figure[3]in large red symbols. The 
analysis of the stacked spectra (red symbols) gives the 
line widths for all LOS within a radial annulus in a lu¬ 
minosity-weighted sense, thus also including faint LOS. 

We find that the line widths of the stacked spectra in¬ 
cluding all LOS are typically larger (by 2 -15%) than the 
values obtained when stacking individual high SNR LOS 

10 This problem does not exist when stacking all LOS. In that 
case we do not perform a preselection, we simply stack all LOS 
falling within the corresponding radial bin. 


only. A possible explanation is that the faint LOS have 
systematically wider line widths than the high SNR LOS 
and when included in the stacks, they widen the result¬ 
ing spectral line width. Since high SNR LOS are less 
than 20% of the total LOS in the CARMA data (which 
account for ~ 60% of the interferometric flux), we can¬ 
not rely on these measurements as being representative 
of the underlying weaker emission. 

4.2. Comparison of FWHM from Different Instruments 

A visual inspection of Figures [3] and [4] already in¬ 
dicates that the CARMA data give the smallest line 
widths of all data sets (in those areas not affected by 
beam smearing). Figurep shows the ratio of the FWHM 
line widths measured by the single-dish telescopes: HE¬ 
RACLES (diamonds) and Nobeyama (squares) divided 
by the line widths measured by CARMA for our two 
galaxies: NGC4736 (Fig. [5] left column) and NGC5055 
( Ki y.tn r ight column). The figure is divided in two. In 
Fig]5(a)| we show the result s for the stacked spectra in¬ 
cluding all LOS. In Fig. 5(b) we show the results involving 
high SNR LOS: on top the stacked spectra, and on the 
bottom the median values of the individual LOS. For 
comparison, we also show the ratio of the line widths 
of the THINGS Hi data (triangles) with the CARMA 
CO (1 0) data. Open symbols show spectra fitted by a 
double - horn profile (typically at small radii where beam 
smearing is dominant), and solid symbols correspond 
to Gaussian fits (typically in the galactic disk]]]] The 
shaded regions mark the galactocentric distance where 
beam smearing is larger than 30% for NGC4736 and 
NGC 5055, respectively. Beam smearing tends to repre¬ 
sent the large - scale gas kinematics instead of the intrin¬ 
sic line profile. As this affects all data sets equally (after 
convolution to a common resolution), it drives any ratios 
of line widths toward unity, thus concealing any intrin¬ 
sic differences. This is obviously the case at small radii 
(inside the shaded regions). At larger radii, where beam 
smearing is not important, the CO line widths measured 
from single - dish data sets (but also for the interferomet¬ 
ric Hi data) are larger than the line widths measured 
from the interferometric CARMA data. 

In Table [ 2 ] we list the mean values for the ratios of 
the single - dish line widths from HERACLES and Nobe- 
maya to the interferometric line widths from CARMA, 
for data only mildly affected by beam smearing, i.e., 
outside 0.15 R 25 for NGC 4736 and outside 0.2 R 25 for 
NGC 5055. The main finding of this paper is that the 
average ratio of single - dish to interferometer line width 
is 1.4 ± 0.2 when looking at the CO (1 - 0) transition and 
taking into account only the high SNR LOS. This value 
is in agreement with the ratio found when taking into 
account all LOS (1.5±0.1). Thus single-dish observa¬ 
tions typically trace molecular gas with larger line widths 
than interferometric observations. This points to the ex¬ 
istence of a molecular gas component that is missed by 
the interferometer and which has larger line widths. 

4.3. Tests 


11 In Figure|5(b)| (bottom) we plot only filled symbols. There is 
no information about line profiles in this plot, since we are showing 
the median of individual measurements. 
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Figure 3. The measured FWHM as a function of galactocentric distance for the four different instruments: CARMA (top left), HERACLES 
(top right), THINGS (bottom left), and Nobeyama (bottom right). The small grey symbols show individual LOS with peak SNR> 5. The 
large red symbols show the values obtained when stacking all LOS within a given tilted ring bin, the medium blue symbols show the results 
when stacking spectra with SNR> 5 only, and the black small symbols show the median values for the grey points, i.e., the median values 
of the individual high SNR LOS measurements. Inside the shaded area beam smearing contributes by more than 30% to the measured line 
width. 
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Figure 4. Stacked spectra by galactocentric distance in bins of 15” width. The spectra for the 4 data sets are shown together, therefore 
fluxes have been rescaled to peak intensities of unity. The black line (labeled ‘C’) is the CARMA data, the red dashed line (‘H’) is the 
HERACLES data, and the green dash-dotted line (‘N’) is the Nobeyama data. For comparison we also show the THINGS Hi data as the 
black dashed line (labeled ‘T’). 


4.3.1. Radial bin filling factor 

As noted before, in the case of NGC 5055 we face in- 
completene ss for the outermost radial bins in this study 
(see Figure|5(a)|). We test how representative the mea¬ 
sured values are, taking into account the incompleteness 
in these points. To do this we take the last complete 
radial bin in NGC 5055, which goes from 45-60" (cen¬ 
tered at ~O.i 5 R. 25 ). We then take a grid of filling fac¬ 
tors going from 90 -10% (i.e., how many individual LOS 
within the radial bin are taken into account). We se¬ 
lect the corresponding percentage of random individual 
LOS inside the radial bin and we stack them. Finally 
we measure the FWHM from the resulting stacked spec¬ 
trum. We perform this same procedure 1000 times, in 
order to see how much the measurement varies by select¬ 
ing random individual LOS. At the end we calculate a 
mean value of the FWHM computed 1000 times for each 
filling factor value and the dispersion among individual 
measurements. 


We compare the “true” FWHM (i.e., the FWHM ob¬ 
tained when stacking 100% of the LOS inside the radial 
bin) to the mean FWHM obtained for each filling fac¬ 
tor. For the CARMA data the variation is at most 2% of 
the original value. For the HERACLES, Nobeyama, and 
THINGS data the variations are at most 2%, 3%, and 
6%, respectively. For these three datasets, the variations 
rise to 5%, 12%, and 4%, respectively in the case where 
the filling factor is of only 10%. Looking at each data set 
individually we note, as expected, that the dispersion in 
the measurements increases with decreasing filling fac¬ 
tor, going from ~ 1km s _1 for the 90% filling factor, up 
to ~4km s' 1 for the 10% filling factor. In all cases, the 
dispersion is of the same order as the errors calculated 
for the stacked spectra (~2-4km s -1 ). 

4.4. FWHM dependance on azimuthal location 

Since the circular velocities have a larger gradient along 
the minor axis (due to projection effects), we also test 
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Figure 5. Ratio of the FWHM computed for the different instruments over the FWHM computed from the CARMA data for NGC 4736 
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10 


Caldu-Primo et al. 


whether the angular distance to this axis affects our 
FWHM measurements. For each galaxy we mask out 
a wedge of X degrees around the minor axis, where X 
ranges from 5 to 40 degrees in steps of 5 degrees. We 
then stack the remaining LOS and compare the FWF1M 
to the value obtained when using no mask. For both 
galaxies and for all angles the differences in FWHM mea¬ 
surements are less than a channel width (<2.5km s -1 ) 
which is 5 -10% of typical FWHM measurements, and 
therefore insignificant. 

4.4.1. Pointing uncertainties 

At this stage we test whether the single - dish pointing 
uncertainties could be responsible for the excess in the 
FWHM measurements. We convolve each single-dish 
cube with a Gaussian of width equal to the correspond¬ 
ing pointing uncertainty^] We compare the resulting 
FWHM to the FWHM measured from the original cubes. 
In all cases the differences between both measurements 
are less than 2%. The pointing uncertainties are thus 
insignificant in the determination of FWHM values. 


4.4.2. Filtering of extended emission 

In this subsection we simulate what the interferometer 
would detect by using the task simobserve in the Com- 
mon Astronomy Softwa re Applications package (CASA) 
(McMullin et al.| 2007). As our galaxy template, we 
use the HERACLEb single-dish cube, which contains 
all emission (clumpy and/or diffuse). We simulate 
what CARMA’sE-configuration UV-coverage would 
recover. We here choose CARMA’s most compact config¬ 
uration since this configuration recovers the largest emis¬ 
sion scales. After constructing the interferometric cube 
using simobserve , we analyze the simulated interferomet¬ 
ric data in the same way as we did with the real data. 
Finally we compare the FWHM measured in the sim¬ 
ulated interferometric observations to the HERACLES 
observations. We find that the HERACLES values are 
larger by a similar amount (~ 50%) to the values mea¬ 
sured in the simulated interferometric data. This further 
supports the idea that the information filtered by the 
interferometer results in narrower line widths, i.e. the ex¬ 
tended emission filtered out by the interferometer has 
also larger velocity dispersion. 

Finally, we caution that CARMA and Nobeyama ob¬ 
served the CO (1 -0) transition whereas HERACLES ob¬ 
served CO (2-1). The mean ratio of the Nobeyama 
to CARMA line widths is ~ 25% larger than the cor¬ 
responding mean ratio of HERACLES to CARMA line 
widths. This difference might be due to different excita¬ 
tion temperatures within the clumpy or diffuse gas or due 
to different optical depth effects. With the available data 
we cannot differentiate between those two mechanisms. 


4.5. Flux Comparison 

We compare the fluxes recovered by the single-dish 
and by the interferometric data sets by means of inte¬ 
grated intensity maps (i.e., zeroth-moment maps). To 
construct these maps we first create 3D - masks that we 

12 The IRAM 30m telescope has a pointing accuracy of ~ 2” and 
the Nobeyama 45m telescope as a pointing accuracy of < 7”. 


Table 2 

Linewidth Ratios 



HERACLES rn (0 n / 

carma 0 “ (2 - o ;> 

Nobeyama co ^ 1 _ 0 ^/ 
carma^q (q 0 ) 

Stacked, all LOS a 

1.2 it 0.1 

1.5 it 0.1 

Stacked, SNR > 5 b 

1.3 ± 0.1 

1.4 ±0.2 

Median, SNR > 5 c 

1.3 ± 0.2 

1.3 ± 0.2 


a Mean of the values measured when stacking all LOS. 
b Mean of the values measured when stacking individual LOS with 
SNR >5. 

c Median of the values measured for the individual LOS with SNR 
>5. 

use to blank the noise dominated regions in the cubes. 
We construct these 3D -masks by finding regions with 
peak SNR larger than 5cr in at least two consecutive 
channels. We then expand these masks to masks con¬ 
structed in the same way, but using a SNR cut of 2cr 
instead. Finally, we expand the mask by half a beam 
size and include adjacent velocity channels to capture 
all emission. We construct two masks for each galaxy; 
one for the interferometric CARMA data set and an¬ 
other one for the single-dish data sets. For the sin¬ 
gle - dish data sets we use the HERACLES cubes to con¬ 
struct the masks, as it has higher SNR as compared to 
the Nobeyama data sets. The final mask is applied to 
each cube and the integrated intensity is calculated. 

We compare the fluxes in three ways: A) We com¬ 
pute the integrated intensity map for each data set in¬ 
dependently (each data set with its own 3D-mask). We 
then measure the flux inside the region limited by the 
CARMA sensitive region (indicated in Figure]!]). B) We 
compute the integrated intensity map for each data cube, 
but use the CARMA 3D - mask in all cases. C) We com¬ 
pute the integrated intensity map for each data cube, 
but use the HERACLES 3D mask in all cases. We then 
measure the flux inside the CARMA - sensitive region. 
The uncertainties we state in Table]!] only take into ac¬ 
count the canonical 10% flux calibration uncertainty for 
each measurement, which is likely a lower limit to the 
actual uncertainties. 

The percentage of flux recovered by the interferometer, 
compared to the flux measured by the single-dish tele¬ 
scope, for the three different approaches are presented 
in Tabled Method A tells us how much flux is recov¬ 
ered by each instrument, even if the flux is originating 
in not exactly the same position, but in the area where 
CARMA has sufficient sensitivity to pick up emission. 
Method B is constrained to look for emission where there 
is CARMA emission, whereas Method C is constrained 
to look for emission where there is HERACLES emission. 
These results show that even with Method B, where we 
measure flux in precisely the interferometric 3D - mask, 
the interferometer recovers ~ 74-81% (depending on the 
galaxy) of the flux recovered by the single-dish. Using 
Method A, which looks for flux in the whole region were 
CARMA has good sensitivity, we find that in the case 
of NGC4736 the flux recovery by the interferometer is 
~ 52% only, while for NGC 5055 it is much higher (92%). 
With Method C, which by construction looks for emis¬ 
sion using the single - dish 3D - mask, the interferometer 
recovers ~ 35% of the flux in the case of NGC 4736 and 
~ 92% in the case of NGC 5055. In the case of NGC 5055 
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Table 3 

Flux Recovery by the Interferometer 3 . 



Method A 

Method B 

Method C 

NGC 4736 

(52 ±14)% 

(74 ±14)% 

(35 ±14)% 

NGC 5055 

(92 ± 14)% 

(81 ±14)% 

(92 ± 14)% 


3 The percentages show the fraction of flux recovered by 
the interferometer (CARMA) as compared to the single-dish 
(Nobeyama).The uncertainties are computed using the typical 
10% flux calibration uncertainty for each measurement. 

the flux recovery is approximately the same, regardless of 
the mask used for the computation of the moment maps. 
In the case of NGC4736 the flux recovery varies from 
~35-74%, depending on the choice of mask. 

Regarding the interferometric flux measurements, we 
recall that during the image deconvolution process, the 
flux recovery will also depend on the noise properties 


(e.g., Heifer et al. 


2002). The cleaning procedure recov¬ 


ers only the flux that is above a few times the noise rms. 
Therefore, the flux that is not being recovered by the 
interferometer due to the noise properties is a low ampli¬ 
tude component. 

Finally, as a consistency check, we compare the flux 
from the Nobeyama 45m cubes to the flux measured 
by the NRAO 12m telescope, observed as part of BIMA 


SONG (Heifer et al. 


2003). For NGC4736 we find that 


the flux measured in Doth observations agrees within 8%, 
whereas for NGC5055 the agreement is within 15%. In 
the case of NGC4736 the interferometer is clearly not 
recovering all the flux measured by the single dish. In 
NGC 5055 the single-dish and interferometric flux mea¬ 
surements agree within their uncertainties. 

5. DISCUSSION AND SUMMARY 

We present a comparison between FWHM line widths 
obtained from interferometric CARMA CO (1-0) data 
to single-dish CO (2- 1) data from the IRAM 30m and 
to single-dish CO (1 0) Nobeyama 45m telescopes in 
two nearby spiral galaxies: NGC 4736 and NGC 5055. 
After convolution of the data sets to a common spa¬ 
tial and spectral resolution, we measure single - dish line 
widths that are ~(40±20)% larger than interferomet¬ 
ric ones (when looking at the common CO (1-0) transi¬ 
tion) and when taking into account the high - SNR - only 
LOS. If we take into account all LOS for the stacking, 
we get single-dish FWHM that are ~(50± 10)% larger 
than the interferometric ones. If the interpretation of a 
diffuse molecular gas component is correct, we would ex¬ 
pect the highest SNR LOS to be dominated by the emis¬ 
sion coming from GMCs. The diffuse emission would be 
more prominent in lower SNR LOS. Furthermore, the 
flux comparisons we perform point to a flux recovery by 
the interferometer as low as ~ 52 - 74% of the flux mea¬ 
sured by the single - dish telescope over the same region 
in the case of NGC 4736. This result is consistent to what 


Pety et al. (2013) found in the case of M51. In the case 
of NGC 5055 the interferometric flux recovery is higher 
~81 -92%, and taking into account the uncertainties in 
the flux determination, we cannot rule out a full flux re¬ 
covery by the interferometer for this galaxy. In summary, 
we measure large line widths and larger flux recoveries 
from the single-dish observations as compared to the 
interferometric ones. The molecular gas ‘unseen’ by the 


interferometer could have two possible configurations. 

One possibility is that this gas is diffuse, i.e., with lower 
average densities as compared to the gas present inside 
GMCs. There is evidence for diffuse molecul ar gas i n 
our Galaxy observed through a bsorption lines (Snow & 
McCall||2006| |Burgli et al.l|2010| 

densities of the diff 


and references therein). 
fuse gas are found to be 


The column densities of t 
< 30 times smaller than those in molecular clouds and 
with typical temperatures d etermined from diffuse H 2 to 
be > 30K ( Liszt et al.|2010 ). This diffuse gas would need 
to extend over volumes larger than the maximal scales 
traced by the interferometer, and thus be filtered out. 
The other possibility is that the gas is inside small molec¬ 
ular clouds which are separated among each other by less 
than a synthesized beam size. This configuration would 
appear as a large regular structure to the interferometer, 
and therefore would be filtered out. The physical proper¬ 
ties of these clouds (temperature and density) may be (or 
not) comparable to the values measured inside GMCs. In 
this case larger velocity dispersion would be the result of 
cloud to-cloud motions. With the available data we 
cannot distinguish between these two possible scenarios. 
Detailed analysis of CO isotopes might help discrimi¬ 
nate between diffuse and dense gas (Liszt et al. 2010), 
as would be higher spatial resolution observations. 

In Galactic studies, the existence of a diffuse molecular 
gas is well accepted (e.g., |Liszt & Pety|2012||Sheffer et al. 
20081. However, these studies are hampered by uncertain 
distances and projection effects that make it difficult to 
quantify the amount of gas present in this pervasive com¬ 
ponent. Studies of nearby galaxies offer an ‘outside view’ 
that do not face these difficulties, though measurements 
are still sparse. Firmly establishing the existence of such 
a spatially extended and high - velocity dispersion molec¬ 
ular gas phase would have important implications. The 
velocity dispersion measured from molecular gas emis¬ 
sion (thought to arise from GMCs only) has been broadly 
used fo r molecular gas scaling relations, or Larson’s laws 
(Larson 1981), among other things, to study star forma¬ 
tion. The presence of two different molecular gas con¬ 
figurations could imply that these scaling relations may 
not be constant throughout the galaxy (depending on 
the dominant molecular gas configuration). Specifically 
in the extragalactic context, where different instruments 
(single dish and interferometer) yield different velocity 
dispersion measurements, studies of a larger sample are 
needed in order to quantify how the scaling relations (for 
example velocity dispersion - molecular cloud size) would 
be affected. 

As a final comment, a complimentary observational ap¬ 
proach to address the same question is to stu dy edge - 
on galaxies, as was done by Garcia-Burillo et al. (1992). 
They detect molecular gas up to ~ 1 kpc oft th± plane 
of the galaxy. They interpret this as evidence of a thick 
molecular disk in NGC 891, which is also in agreement 
with a diffuse molecular gas component, since higher ve¬ 
locity dispersions in principle imply a larger scale height 
of the molecular gas disk. Future observations with 
ALMA, using both interferometric observations as well 
as total power measurements, will be crucial to shed more 
light on the nature of this ‘diffuse’ molecular gas phase 
in galaxies. 
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